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Abstract 

We study the structure-dependent contributions to the radiative baryonic B decays of B — > BB'7 
in the standard model. We show that the decay branching ratios of Br(B — > BB'7) are O(10~ 7 ), 
which are larger than the estimated values of O(10~ 9 ) induced from inner bremsstrahlung effects of 
the corresponding two-body modes. In particular, we find that Br(B~ — > Apj) is around 1 x 10~ 6 , 
which is close to the pole model estimation but smaller than the experimental measurement from 
BELLE. 

PACS numbers: 
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The radiative baryonic B decays of B — > BB'7 are of interest since they are three- 
body decays with two spin- 1/2 baryons (B and B') and one spin-1 photon in the final 
states. The rich spin structures allow us to explore various interesting observables such as 
triple momentum correlations to investigate CP or T violation . Moreover, since these 
radiative decays could dominantly arise from the short- distance electromagnetic penguin 
transition of b — > S7 yl which has been utilized to place significant constraints on physics 
beyond the Standard Model (SM) |4J,|5(, they then appear to be the potentially applicable 
probes to new physics. 

There are two sources to produce radiative baryonic B decays. One is the inner 
bremsstrahlung (IB) effect, in which the radiative baryonic B decays of B — * BB'7 are 
from their two-body decay counterparts of B — ► BB' via the supplementary emitting pho- 
ton attaching to one of the final baryonic states. Clearly, the radiative decay rates due to 
the IB contributions are suppressed by a em comparing with their counterparts. According 
to the existing upper bounds of B — > BB', given by 

Br(B° -> pp) < 2.7 x 1(T 7 (BABAR) , 
Br(B° AA) < 7.9 x KT 7 (BELLE), 
Br(B~ -> Ap) < 4.6 x KT 7 (BELLE) , 



(1) 



one finds that 



Br(B -> BB'7) /B < O(10" 



(2) 



Unfortunately, the above branching ratios are far from the present accessibility at the B fac- 
tories of BABAR and BELLE. However, the other source, which is the structure-dependent 
(SD), is expected to enhance the decays of Br(B — > BB'7), such as B — > Ajry arising from 
— >• S7 . With the large branching ratio of b — > S7 

expect that Br[B~ BB'7) could be as large as Br(B~ - 
concentrate on the SD contributions to Br(B — > BB'7). 
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12| in the range of 10 4 we 



BB'). In this report, we shall 



To start our study, we must tackle the cumbersome transition matrix elements in B — > 

14L I15| in recent years, 



BB'. As more and more experimental data on three-body decays 



the theoretical progresses are improved to resolve the transition matrix element problems. 



One interesting approach is to use the pole model 
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and another one is to rely on the QCD counting rules [18 



irough the intermediated particles 
^a by relating the — 



FIG. 1: Diagram for B — > Ay>7 

matrix elements with three form factors and fitting with experimental data. In Ref. j^]], 
Cheng and Yang have worked out the radiative baryonic B decays based on the pole model. 
In this paper, we handle the transition matrix elements according to the QCD counting 
rules. 

We begin with the decay of B~ — > Ajry. As depicted in Fig. in the SM the relevant 
Hamiltonian due to the SD contribution for B~ — > Apj is 



n SD = -^V tb V*c e / f 7 , (3) 



with the tensor operator 



7 = ^"1^^(1 + 75)6, (4) 



where V t bV t * and cy are the CKM matrix elements and Wilson coefficient, respectively, and 
the decay amplitude is found to be 

A(B- - Am) = ^=V t lV tb ^- 2 2c e / f 

mge"<Ap|s 7M (l - j s )b\B~) - 2m b p B ■ e(Ap\s(l + j 6 )b\B-)\, (5) 



3 



where we have used the condition m& ^> m s such that the terms relating to m s are neglected. 
We note that Eq. (jSJ) is still gauge invariant. 

In order to solve the encountered transition matrix elements in Eq. we write the 

most general form 

(Ap^^B') = m(p A )[ai7 M 7 5 + a 2 p^ 5 + a 3 (p P - P\)lb\v{pp), 
(Ap 1 57^756 \B~) = iu(p A )[c 1 'j IM + c 2 ia l , u p u + c 3 (pp + p A )]v(p p ), (6) 

where p = Pb — Pa — Pp and cii(ci) (i = 1, 3) are form factors. 

To find out the coefficients Oj(cj) in Eq. ©, we invoke the work of Chua, Hou and Tsai in 
Ref. j^J. In their analysis, three form factors Fa , Fp and Fy are used to describe B — > BB' 
transitions based on the QCD counting rules (l8[, that require the form factors to behave 
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as inverse powers of t = (p& + Pb>) 2 - The detail discussions can be referred to Refs. 
In this paper, we shall follow their approach. The representations of the matrix elements 
for the B~ — > pp transition are given by jscj 

(pp\u(l ± j 5 )b\B~) = iu{p p )[{F A ^ 75 ± F v -f) + (F p75 ± F s )]v(p p ), (7) 



with a derived relation F$ = Fp. In terms of the approach of |19Ll20|. those of the B — > Ap 
transition are given by 

(Ap\s(l ± 7fe)6|B-> = m(p A )[(Ff ^75 ± F? f) + {F** 75 ± Ff)]v{pp), (8) 
where the form factors related to those of B~ — > pp in Eq. J7J) are shown as 

The three form factors Fa, Fy and Fp can be simply presented as 19 



f*.v = % = % do) 

where Cj (i = A, V, P) are new parametrized form factors, which are taking to be real. 

From the relation p /i (Ap|s7 At (l — 75)6|.B _ ) = m&(Ap|s(l — 7s)6|i?~) in the heavy b quark 
limit, the parameters dj(cj) in Eq. (JBJ) are associated with the scalar and pseudo-scalar 
matrix elements defined in Eq. (JBJ). As a result, we get that 



7-iAp m bFp P Ap m bFp P , , 

ai = m b F/ , a 3 = — , c x = m b F v F , c 3 = — ■ -. (11) 



P ■ (Pp ~ Pa) ' P ■ (Pp + Pa) 
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The amplitude in Eq. @ then becomes 
A(B~ 



Gi 
V2 



fTw + FiH ipp Pa) " 



V-ijPp-Px) V M P V ■ (Pp + Pa) 



v (Pp 



-2m hPB ■ e u(p A ) Ff fa + F* p l5 + F^ j + F 1 ? v(p p 



Ap 



? Ap 



^Ap 



(12) 



with three unknown form factors F^ p , Fy p and Fp p . We note that the terms corresponding 
to a 2 disappear due to the fact of e ■ p = 0. Even though c 2 can only be determined by 
experimental data, according to QCD counting rules, c 2 needs an additional 1/t than c\ to 
flip the helicity, so that it is guaranteed to give a small contribution and can be neglected. 

After summing over the photon polarizations and baryon spins, from Eq. ()12j) . the decay 
rate of T is given by the integration of 



dT 



where 



(2tt) 3 4M|£2 



\C t 



V\F Ap \ 2 + A\F Ap \ 2 + P\F Ap \ 2 + I VP Re(F Ap F Ap * 



+I AP Re(F Ap F Ap * 



dm\ p dm 2 pi 



(13) 



rriAp = PA+Pp 



m 



pi 



G F 

Pp + p-,, C t = —7=V tb V t 



V2 



-2c 



eff 



V(A) = p A ■ p(E p E y - pp ■ p) + E 7 (Eap p ■ p ± E y m A m p ) 



£ 7 (E A + Ep) {m A mp - p A ■ Pp) (m\ + m 2 p + 2p A ■ p P ) (m A mp - p A ■ Pp) 



PA-P + Pp-P 
E 7 (E A - Ep)(m A mp - p A ■ p p ) 



PA-p-Pp-p 



2(p A ■ p + Pp ■ p) 2 
(m\ + m? - 2p A ■ Pp) (m A mp + p A ■ p p ) 
2 (pa -p-Pp-p) 2 



-PA ■ Pp 



Wp(ap) = 2EpE 1 rriA - Pp • pm A ± E A Ej(mA - m p ) ± mpp A ■ p 

E^(E P ± E A )(m A + m p )pA ■ p - E 2 (m A - m p )(p A ■ p p ± m A mp) 

Pa ■ p ± Pp ■ p 



(14) 



It is important to note that, since the penguin-induced radiative B decays are associated 
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with axial- vector currents shown in Eq. (jSJ), we have used 21] 

E*a a _ , k^n v + k u n 

A=l,2 



k ■ n 



k^k u 
(k ■ n)'- 



(15) 



where n = (1, 0, 0, 0), to sum over the photon polarizations instead of the direct replacement 



A=l,2 C M c u 



~9fiu which is valid in the QED-like theory due to the Ward identity. 



For the numerical analysis of the branching ratios, we take the effective Wilson coefficient 



= —0.314 |22j, the running quark mass nib = 4.88 GeV and CKM matrix elements 
V t bV t * = —0.0402. Even though there are no theoretical calculations to the unknown Ca, Cy 
and Cp. By virtue of the approach of Ref. j^J, these form factors are related to the present 
experimental data, such as Br(B~ — > ppn~), Br(B° — > ppK°), Br(B~ — ► ppK~) (l5j and 



Br{B~ 



AAK~ 



characterized by an emitted pseudoscalar meson. For a reliable % 2 
fitting, we need 2 degrees of freedom (DOF) by ignoring the Cp term since its contribution 
is always associated with one more 1/t over Ca and Cy ones, as seen in Eq. (jlOj) . We will 
take a consistent check in the next paragraph to this simplification. To illustrate our results, 
we fix the color number Nq = 3 and weak phase 7 = 54.8°. The input experimental data 
and numerical values are summarized in Table HJ 

TABLE I: Fits of (C A ,C V ) in units of GeV 4 . 



Input 


experimental data 


Fit result 


best fit (with 1 a error ) 


Br{B~ -> ppTT") [15] 
Br(B° -» ppK°) [15J 
Br{B~ -^ppK~) [15] 
Br{B~ -> KAK-) [23] 


3.06 ±0.82 
1.88 ±0.80 
5.66 ±0.91 
2.91 ±0.98 


C A 
Cy 
X 2 /DOF 


-68.3 ±5.1 
35.1 ±9.0 
1.85 



Using the fitted values of C A and C Vl we find BriB' Apj) = (0.92 ± 0.20) x 10" 
which is larger than its two-body decay partner as expected and it is close to the result of 
1.2 x 10~ 6 in the pole model 



than (2.16^;1 ± °- 20 ) x 10 ' 



However, our predicted value on B~ — > A/ry is smaller 
2J] measured by BELLE. If we put this new observed value 
into our fitting, we can further include Cp ignored previously. The fitted values are Ca = 
-73.3±9.1 GeV\ Cy = 43.7±12.1 GeV 4 and C P = 134.3±327.0 GeV 7 with x 2 /DOF = 3.65 
which is about two times bigger than previous one. Clearly, it presents an inferior fitting 
with small Ca,v changes. When putting back these three fitted values to Br{B~ — > Ajry) for 
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2.5-10 




FIG. 2: dBr(B — > Ajry) / dm\p vs. mxp- The solid line stands for the input values of (Ca, Cy) = 
(-68.3,35.1) while the dash line stands for those of (C A ,C V ,C P ) = (-73.3,43.7, 134.3). 

a consistency check, we get (1.16 ± 0.31) x 10~ 6 regardless of inputting larger experimental 
value, which explains the large value of x 2 /DOF. The insensitivity of Cp on the decay 
branching ratio justifies our early simplification of ignoring its contribution beside the l/t 
argument . 

In Ref. Q, it was SUgg e s ted that the .educed e Mgy ^ ean mak e the branchy 
ratios of three-body decays as significant as their counterparts of two-body modes or even 
larger, and one of the signatures would be baryon pair threshold effect |l|, |20j. In Fig. |2l 
from Eq. (|T3*|) we show the differential branching ratio of dBr(B~ — > Ap r y)/dm^p vs. m\p 
representing the threshold enhancement around the invariant mass m\p = 2.05 GeV, which 



is consistent with Fig. 2 in Ref. 



24j of the BELLE result. Around the threshold, the 



baryon pair contains half of the B meson energy while the phone emitting back to back 
to the baryon pai r with another half of energy which explains the peak at E n ~ 2 GeV in 



Fig. 3 of Ref. 
method works 



2M- Such mechanism is similar to the two-body decays so that factorization 
1( even in the three-body decays. 
To discuss other radiative baryonic B~ decays, we give form factors by relating them to 
Fy,A,p m the B~ — > pp transition similar to the case of B~ — * Ajry as follows: 



B- 



llFy 
10V2 



10V2' A 



9F V 



11F A 



F Z 0p = _t^_ 

10\/2 lOv^' P 3\/2' 
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B — > X 777 : 

llFy 



V 



10 



lb - ' F ^ 



9F y 9F A 



10 



11 



Fp~ n 



Fp 

4 : 



5" 



5 A7 : 



lOv 7 ^ lOvV A 



lOv^ lOvV 



F P 

,*p --j-, 



B- 



>E- 7 : 
10 " 



"10"' FA 



9^V _ ^ s o g - 5Fp 
10 10' p 4 ' 



B' 



9F 



A p3 -S° 



9F V 



5-Fc 



(16) 



10v^ lOv 7 ^' 10^2 10>/2' 4^2' 

To calculate the branching ratio of B — > BB' 7 , we can use the formula in Eq. (|13|) by 
replacing A and p by B and B', respectively. The two sets of predicted values for B — > BB'7 
with and without Cp are shown in Table ITTl respectively. As a comparison, we also list the 
work of the pole model approach by Cheng and Yang in the table. We note that, in 



TABLE II: Decay branching ratios 





Fits 




Branching Ratios 


(Ca, Cy) = 


(Ca, Cy,Cp) — 


Pole model[9j 




(-68.3 ±5.1, 35.0 ±9.0) 


(-73.3 ± 9.1, 43.7 ± 12.1, 134.3 ± 327.0) 




Br(B~ -> Ap 7 ) 


(0.92 ±0.20) x 10~ 6 


[(1.16 ±0.31) x 10~ 6 ] 


1.2 x 10~ 6 


Br(B~ -> E°p7) 


(1.7 ± 1.5) x 10~ 7 


(1.2 ± 1.2) x 10~ 7 


2.9 x 10~ 9 


Br(B- -> S"fi7) 


(3.4 ±2.8) x 10~ 7 


(2.5 ±2.4) x 10~ 7 


5.7 x 10~ 9 


Sr(S" -> H-A7) 


(0.48 ±0.50) x 10" 7 


(0.61 ±0.60) x 10" 7 


2.4 x 10~ 7 




(3.3 ±0.7) x 10~ 7 


(3.7 ±0.9) x 10~ 7 


1.2 x 10" 6 




(1.5 ±0.6) x 10~ 7 


(1.8 ±0.6) x 10~ 7 


6.0 x 10~ 7 



Table ITTl the value in the bracket of the third column for Br(B~ — > Apy) is not a prediction 
but a consistency comparison with the puting-back form factors, since we have used the 
observed value of Br(B~ — > Ap~y) from BELLE. We found that, except for Br(B~ — > Apy), 
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all predicted values are O(10" 7 ). In terms of inverse sign between Ca and Cy, there are 
constructive effects for and F v ^, which are proportional to {Fy — Fa) as shown in Eq. 
(JSJ), whereas destructive effects make other F^ B and Fy B ' in Eq. (J16)) small. Consequently, 
all modes for B~ radiative baryonic decays are suppressed except for Br(B~ — > Afry). We 
remark that such suppressions exist only in the SM-like theories. Thus, these radiative 
baryonic decays are useful modes for testing the new physics. 

As seen in Table |HJ both our results and those of the pole model satisfy the relations 
of Br(B~ Trni) ~ 2Br(B~ -> S°p7) and Br(B~ -> H £-7) ~ 2Br(B~ S~S°7) 
because of the SU(3) symmetry. In the pole model, the decay branching ratios of B~ — > Apj 
and B~ — > H°X~7 are found to be large, around 1.2 x 10~ 6 , since ther are intermediated 
through A b and E b , which correspond to large coupling constants g\ b ^B-p an d fl , sO-^s-s+, 
respectively. However, in our work, the branching ratio of B~ — > Ap7 is about three times 
larger than that of B~ — > S°S~7, which is O(10~ 7 ). Regardless of these differences, both 
two methods are within the experimental data allowed ranges, such as those of 

[Br(B- -> Ap7) + 0.3Br(B- -> £°p 7 )]£ 7>2 . GcV < 3.3 x 1(T 6 , 
[Br(B- -> S°p 7 ) + 0.4Br(B- -> Ap7)]s 7 >2.o Gov < 6.4 x 1(T 6 , 



from CLEO [25| and 5r(B~ -> S p 7 ) < 3.3 x 1(T 6 from BELLE |24|. 

Finally, we relate the B° decays with the corresponding B~ modes in terms of QCD 
counting rules even though there are no experimental data on radiative baryonic B° decays. 
When neglecting the mass and life time differences, we obtain 

Br{B~ -> Ap7) = Br(B° -> Arry), Br(B~ -> £°P7) = Br(B° -> S n7), 
£r(-B _ -> £~n 7 ) = Br(B° -> S + P7), Br(B~ -> S-A7) = £r(£° -> S°A 7 ), 
Br(S~ _> S"S° 7 ) = Br(S° -> S°S°7), Br(S" -> S°S- 7 ) = fir^ -> --£+7) , (17) 

which are also guaranteed by the SU{3) symmetry. From Eq. 1)1 7]). we see that Br(B° — > 
A717) can be as large as Br(B~ — > Ajry). 

In sum, we have shown that the SD contributions to the radiative baryonic decays of 
£> — > BB'7 in the SM are associated with the form factors of Fa, Fy and Fp in the matrix 
elements of the B~ — > pp transition. Most of the predicted values for Br(B — > BB'7) 
are spanning in the order of 10~ 7 , which are larger than the estimated values of O(10~ 9 ) 
due to the IB effects of their two-body counterparts. In particular, we have found that 



Br(B~ Ap7) is (1.16 ± 0.31) x 10~ 6 and (0.92 ± 0.20) x 10" 6 with and without C P , 
respectively, which are consistent with the pole model prediction [l| but smaller than the 
experimental data from BELLE 2^|. More precise measurements are clearly needed. 
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